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Abstract. Wc study the cosmic ray arrival distribution expected from a source 
of neutrons in the galactic center at energies around 1 EeV and compare it with 
the anisotropy detected by AGASA and SUGAR. Besides the point-like signal in 
the source direction produced by the direct neutrons, an extended signal due to 
the protons produced in neutron decays is expected. This associated proton signal 
also leads to an excess in the direction of the spiral arm. For realistic models of the 
regular and random galactic magnetic fields, the resulting anisotropy as a function 
of the energy is obtained. We find that for the anisotropy to become sufficiently 
suppressed below E ~ 10^^'^eV, a significant random magnetic field component is 
required, while on the other hand, this also tends to increase the angular spread of 
the associated proton signal and to reduce the excess in the spiral arm direction. 
The source luminosity required in order that the right ascension anisotropy be 
4% for the AGASA angular exposure corresponds to a prediction for the point- 
like flux from direct neutrons compatible with the flux detected by SUGAR. We 
also analyse the distinguishing features predicted for a large statistics southern 
observatory. 
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1. Introduction 

The distribution of the arrival directions of high energy Cosmic Rays (CRs) provides 
a useful tool to investigate their origin. This distribution is found to be very close 
to isotropic up to the highest energies measured. For energies below the knee 
{Eknee — 3 X 10^^ cV) the measured anisotropy is below 0.1%, showing a smooth 
increase for larger energies and approaching the percent level around 10^^ eV 
This isotropy is understood as due to the diffusive motion of charged particles in 
the galactic magnetic field for energies such that E/Z < 10^^ eV, with Z the CR 
charge. Cosmic rays from galactic sources remain trapped in this field, leading to a 
flux building up at small energies and to the absence of correlations between the arrival 
direction of CRs and the source's direction. At energies above 10^^ eV, and up to the 
highest observed energies, the distribution of the events is compatible with isotropy. 
However, at these energies CRs are no more trapped in the galactic magnetic field, so 
that a correlation between the arrival and the source directions is to be expected. The 
absence of a significant correlation between ultra high energy events and the Galactic 
plane points then to extragalactic sources as the origin of these events. 

In 1999 the AGASA group reported the measurement of the most significant 
anisotropy detected so far [2] . They presented the results of a search for anisotropics 
in the arrival directions of cosmic rays with energy above lO" eV, using data collected 
over 11 years. They found a strong anisotropy of amplitude ~ 4% in a small energy 
range: 10^^'^ — 10^* '^ eV. The phase of the first harmonic in right ascension was 
found to be ai ~ 300°. A two dimensional map showed that this anisotropy can be 
interpreted as due to excesses of events in two regions of size ~ 20° near the galactic 
center (GC) and the Cygnus region with significance 4cr and 3a, respectively. A 
subsequent analysis by the same group, including an independent data set, confirmed 
the initial claim 0. 

Taking into account that the excess reported by AGASA lies close to the galactic 
center, that is better seen from the south, the SUGAR group reanalyzed their data 
The SUGAR experiment, located near Sidney, has a clear view of the GC, but has 
much poorer statistics than AGASA. They set a priori the energy range corresponding 
to the AGASA anisotropy and found a point like excess at 7.5° from the GC, and 6° 
away from the AGASA maximum. Due to the difference in direction and angular 
extension of the AGASA and SUGAR signals it is not yet clear whether they are 
consistent among each other and if they can result from a common origin. 

A possible explanation of the anisotropy, already proposed by Hayashida et al. 
12], is that it is due to neutron primary particles. Neutrons of 10^* eV have a Lorentz 
factor 7 ~ 10®, so that their decay length is about 10 kpc, i.e. similar to the distance 
from the Earth to the galactic center. Hence, neutrons from a source near the galactic 
center can propagate up to the Earth before decaying, and they are undeflected by 
intervening magnetic fields, contrary to the case of charged primaries. The lack of 
significant anisotropics at smaller energies is explained within this scenario as due to 
the shorter decay length of the neutrons, and at larger energies by a cutoff in the 
acceleration energy (or due to a very steep spectrum) of the sources. They suggested 
that neutrons can be produced by accelerated heavy nuclei which disintegrated due 
to interactions with ambient matter or photons surrounding the source. Charged 
particles could stay confined by the magnetic fields present in the acceleration region, 
while neutrons would escape easily. 

Medina- Tanco and Watson _6_ proposed that a more likely way of generating a 
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high energy neutron flux is through high energy protons via interactions with ambient 
protons or with infra red photons. They found that the neighborhood of the Galaxy's 
central supermassive black hole could give the environment needed for these reactions 
to happen frequently enough. Nagataki and Takahashi showed that the proton-proton 
process was the most efficient one, and also studied the detectability of the 7 and 
emission associated to the production of neutrons in the GC . 

The possibility that the GC is a source of high energy protons was studied by 
Levinson and Boldt |H]. They showed that, under certain conditions, a compact black 
hole dynamo associated with the Sgr A* source may accelerate protons up to energies 
~ 10^® eV, with a total luminosity of ~ 10"^^ erg/s. 

In this work, we study in detail the CR arrival direction distribution expected 
from a source of neutrons in the GC. The fact that the decay length of 10^* eV neutrons 
is similar to the distance to the proposed source implies that a sizable fraction of the 
neutrons will actually decay to protons before reaching the Earth, and the trajectories 
of these secondary protons will be bent by the galactic magnetic fields (GMF). These 
protons would arrive to the Earth from some preferred directions in the sky: those 
produced near the Earth would come from directions close to the GC one, while those 
produced in the inner Galaxy region would arrive preferentially from directions close 
to the spiral arm, as their trajectories wind around the regular magnetic field lines. 
Initial work in this direction was performed by Medina- Tanco and Watson |S] and 
here we extend it to take into account different allowed magnetic field strengths, we 
map the arrival direction distribution for different energies and take into account the 
angular exposure of experiments at different latitudes to study several features of the 
observations such as the energy dependence of the anisotropics, the compatibility of 
the observations of AGASA and SUGAR both regarding the angular extent of the 
observed features and the required source luminosity. 

Other suggested explanations for the observed anisotropics include a diffusive flux 
of charged particles from the GC direction [HI , the drift of heavy nuclei in the galactic 
magnetic field ^01 and a source of protons in the galactic center [111 IT^ . 

2. Cosmic Ray propagation in the Galaxy 

The trajectories of charged CRs will be determined by the galactic magnetic field, 
which has both a large scale regular component and a random component, both with 
strengths of a few fiG. The magnetic lines of the regular field follow the spiral structure, 
reversing direction between consecutive arms. According to the preferred model for 
our galaxy, it is symmetric with respect to the galactic plane. The strength of the 
field can be modeled as I14j 



where {p, 9, z) are cylindrical coordinates with origin in the galactic center, — 10.55 
kpc is the value of p for which the field is maximum in our spiral arm, (3=1/ tan p = 
—5.67, where p = —10° is the pitch angle and tq = 8.5 kpc is the distance from the 
Sun to the GC. The overall intensity is taken as So = 3 /iG, corresponding to a local 
strength of the regular component close to 2 pG. The tanh term smoothes the 1/p 
behavior of the field near the center within a core radius taken as pi = 2 kpc. The 
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vertical scale length associated to the disk field is taken as zi = 0.3 kpc, while the 
halo one as Z2 = 4 kpc. 

The radial and azimuthal components are given by Bp = Bgp sin p and 
Bg = Bap COS p respectively, while the strength of the field in the z direction is 
assumed to vanish. 

Besides the regular magnetic field structure, a random component due to 
turbulence in the interstellar plasma is known to exist, with the largest eddies having 
a scale of Lmax — 100 pc. A usual assumption is that the spectrum of the field 
inhomogeneities is the same as that of the gas density, that is close to a Kolmogorov 
spectrum. Hence, the turbulent magnetic field can be described by a random field with 
a power law spectrum such that the magnetic energy density satisfies dE/dk oc k^^/^ . 
It can then be written as a sum of Fourier modes 

B-^{S) - I 4^BS)e^^^-'^+^^^^\ (2) 



(27r) 

where the phases (j)i{k) are random and the amplitude of the modes is given by 

R2^L^ 2 2 . 5/3 {2T:/Lmaxf^ 

^ {l^) - o^rmsl^ 1 _ iT IT ^2/3' ^'^> 

\^ min I ^ max } 

for 27: / L,-aax ^ k < 27r/L,„i„, and zero otherwise 16 . Brms fixes the root mean 
square amplitude of the random component. In the numerical calculations the integral 
in Equation jSJl was replaced by a sum over 200 modes. Each B{ki) had a random 
amplitude extracted from a Gaussian distribution with zero mean and root mean 
square value given by the square root of Equation Q and the direction was chosen 
randomly in the plane perpendicular to fci, ensuring that V • B{x) — 0. Since particles 
reaching the Earth at the energies under study never departed from the galactic 
plane by more than a few hundred pc, the detailed vertical profile of the random 
component is not relevant for the results, and hence a constant value for it was adopted 
throughout. Moreover, the random field was also taken as having a constant rms 
strength in the radial direction. 

The effect of these field components on the motion of EeV protons can be 
estimated as follows. The gyroradius of a high energy CR {v k, c) with charge 
Ze and energy E traveling through an uniform magnetic field of strength B is 
i?kpc — ^-^E-EcY / ZBpQ. Hence, a proton with E ^ \ EeV in a B '--^ 3 /iG field 
has a gyroradius R ^ 0.4 kpc. As the regular GMF is uniform over scales of the order 
of a few kpc and the thickness of the galactic disc is 1 kpc, the trajectory of an EeV 
proton is expected to wind around the magnetic field lines. 

The turbulent component adds random defiections to the propagation. The rms 
defiection of a CR with energy E propagating a distance L through a turbulent 
magnetic field of strength B^ms and maximum inhomogeneity scale Lmax (for L ^ 
L„iax) is given by ^ 



, EeV Brms I L 



^' s/z 1 ^Gy 10 kpcy 100 pc ^ 

Thus, we expect that for EeV protons traveling in the galactic turbulent field, this 
will induce a significant dispersion in the arrival directions. 
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3. Results 

3.1. Arrival direction distribution 

Neutrons emitted from the GC travel along straight lines until the moment in which 
they decay into protons, and it is only the trajectories of these last which start to 
be bent by the GMF. Thus, in order to obtain the distribution of arrival directions 
at Earth we should follow the trajectories of protons injected in spheres of various 
radius around the GC and with velocities pointing out of the GC, and record the 
direction of those arriving to the Earth. These would have to be weighted by the 
probability of neutron decay at the corresponding radius. However, a more efficient 
way to compute this distribution is to back track antiprotons leaving the Earth in 
all different directions and record the points where their velocities point toward the 
GC. It is clear that if a neutron from the GC were to decay at that point, the proton 
produced would arrive to the Earth. We actually considered a source with finite 
radius, and obtained the distance from the GC (r) and the length of the segment 
(Ar) of the proton trajectories for which the velocity pointed towards the source. The 
corresponding flux of protons from neutron decays that arrive from a given direction 
is then dJ/df2 oc e"''/^'^'^" Ar/7cr„. Adding finally the contribution from neutrons 
which did not decay, the arrival distribution map can be constructed. 

The left panels in Figure ^ show the arrival direction distribution obtained in 
this way when only the regular magnetic field component is included and for different 
values of the energy between 10^^'^ eV and 10^^'^ eV. The most intense flux arrives 
from the GC direction {aac — 266.4° and 6gc = —28.9°) and corresponds to the 
direct neutrons. A thin (nearly horizontal) strip extending from the GC towards the 
north and south galactic poles (i.e. perpendicular to the galactic plane, shown as 
a dashed line) is due to neutrons decaying near the Earth in our spiral arm and in 
the nearest inner arm respectively. The more extended excess appearing for right 
ascensions between 300° and 360° correspond to neutrons that decayed in the inner 
galactic region, with the protons reaching the Earth through the spiral arm. We 
can see that the region of the sky where events are spread shrinks considerably with 
increasing energy, in particular, at 10^^'^ eV almost all the events arrive outside the 
field of vision of AGASA (S > -24°). Notice that the flux intensity of the plots with 
different energy should not be compared at this point, as the energy spectrum of the 
source has not yet been introduced. 

The inclusion of the turbulent magnetic field component leads to a spreading in 
the arrival direction distribution as illustrated in the right panels of Figure^ in which 
Brms = 3 /iG was adopted. The spread in the arrival directions increases with the 
strength of the random component, as indicated by Equation and this is clearly 
seen in the left panels of Figure El which show this distribution for a fixed energy of 
10^^ eV and increasing intensity of the random field {B„ns = 0, 1, 2 and 3 ^G). 

The right ascension density distribution for the Brms = 3 /iG random field maps 
of Figuren]is shown in the left panels of Figure|3| The sharp peak in a ~ 265° is mainly 
due to the direct neutrons. This is relatively more important as the energy is increased, 
due to the increase in the neutron's decay length. The broader excess associated to the 
secondary protons arriving along the spiral arm becomes strongly suppressed above 
10^^ eV due to the increasing difficulty for the regular field to confine the CRs along 
the spiral arms. This suppression is clearly expected to be more important as the 
value of Brms becomes stronger, while on the other hand a larger value of the regular 
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Figure 1. Arrival direction distribution in celestial coordinates for increasing 
values of the CR energy. In the left panels only the regular galactic magnetic field 
is included while in the right panel a random field component with Brms = 3 fiG 
is included. The dashed line represents the galactic plane with the black dot in it 
corresponding to the direct neutron flux from the GC. Darker points correspond 
to greater density in logarithmic scale. 



component Bq should lead to a signal more concentrated in the spiral arm direction 
and surviving up to larger energies. 

In order to compare our results with AGASA's measurements, it is necessary 
to take into account the exposure (lo) of the observatory, which gives the time- 
integrated effective collecting area from each direction of the celestial sphere. The 
arrival direction distribution measured by an observatory should be proportional to 
the actual distribution multiplied by uj. 

The exposure is expected to be homogeneous in right ascension (a) and to have 
a declination (S) dependence given by 

a;((5) cc cos(ao) cos(J) sin(am) + a„i sin(ao) sin(5), (5) 
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Figure 2. The left panels show the arrival direction distribution in celestial 
coordinates for different values of turbulent field strength Brms = 0, 1, 2 and 3 fxG. 
In all plots the energy is E = 1 EeV. The right panels show the corresponding 
right ascension distribution taking into account the AGASA exposure. 



where is given by 



and 



if C > 1 

TT if^<-l 
cos" ^ (^) otherwise 



^ _ cos.{9m) — sin(ao) sin((5) 



cos(ao) cos((5) 

where ao is the detector latitude (in the case of AGASA ao = 35.78°) and 6m is the 
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Figure 3. Right ascension density distribution corresponding to the a random 
field with Brms = 3 fiG for different values of the CR energy. The left panels are 
for the full sky maps, while in the right panels the exposure of AGASA is used. 
The scale is in arbitrary units. 



maximum zenith angle observed. The analysis of AGASA data correspond to events 
with zenith angle < 60°. 

Taking into account the exposure, the right ascension distribution expected to be 
measured at the AGASA latitude looks as shown in the right panels of Figure |21 for 
a random field amplitude Bj-ms = 3 /iG. The main difference with the full sky results 
displayed in the left hand panels is the disappearance of the point-like density peak 
due to the direct neutrons. The right panels of Figurc[21also show this distribution for 
different values of Brms and for E = lO^^cV, for the AGASA exposure. The increasing 
spread of the signal associated to the spiral arm, which is the dominating one in this 
case, is apparent. 

Anisotropics in the CR distribution are often measured performing an harmonic 
analysis of the right ascension distribution of the events This takes advantage 

of the uniform exposure in right ascension of most experiments. The right ascension 
distribution is related to the differential spectrum of CR, dJ/dil{E) through 



F(a) = / dScosSLj(S) — . 
J ail 
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Expanding it as F{a) = Fq + Fi cos(q! — ai) + . . . , the amplitude of the anisotropy 
is defined as the ratio of the dipole to the monopole intensities: A = Fi/Fq. In 
Tabled the phase of the first harmonic (ai) from the analysis of the right ascension 

Table 1. Phases of the first harmonic from the analysis of the right ascension 
distribution as a function of energy and turbulent magnetic field strength, for the 
exposure at the AGASA latitude and at 35° south. 
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distribution are listed as a function of energy and turbulent magnetic field strength, 
taking into account AGASA's exposure, and also for a hypotetical detector in the 
southern hemisphere (at a latitude —35°, similar to the one of AUGER or Sidney). 

3.2. Amplification 

The total flux that we receive from a given source is modified by the presence of a 
magnetic field. This can be described by introducing the amplification, A(E), which 
is the ratio of the flux of particles arriving to the Earth from any direction in the 
presence of the GMF, J{E), and that in the absence of magnetic field, Jq{E), i. e. 
A{E) — J{E)I Jq{E). The amplification is a function of the energy and the mean 
effect turns out to be larger as the energy decreases, as expected since the CRs tend 
to be trapped more efficiently in the spiral arms of the GMF. At larger energies the 
amplification tends to unity, as protons are less affected by the magnetic field and a 
larger fraction of the neutrons arrive to the Earth before decaying. Figure0]shows the 
amplification vs. energy for adopted values of the random magnetic field of Brms ~ 
and 3 /iG. The smaller amplification obtained in the presence of the random magnetic 
field is due to the enhancement of the escape of CRs from the spiral arms produced 
by this component. 

3.3. Anisotropy 

The anisotropy detected by AGASA arises in a narrow range of energies from 
E = 10^^'^eV to E — 10^^ "^eV. It has been suggested that the lack of anisotropics at 
smaller energies could be explained in the model discussed here because of the shorter 
decay length of the neutrons and the deflection of the resulting protons in the galactic 
magnetic field. In order to test quantitatively this hypothesis we have to estimate the 
total anisotropy produced by the GC neutron source in the presence of the background 
of all the other sources of CR in that energy range. 

The fiux of CRs with energy around 1 EeV arriving to the Earth is thought to 
come from galactic sources, such as supernovae, that give the main contribution below 
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Figure 4. Amplification vs. energy of a neutron source located at the GC due 
to the GMF, with (solid line) and without (dotted line) including the turbulent 
component of the GMF {Brms = 3 /iG) in the simulation. 



the ankle {Eankie — 5 x 10^^ eV), as well as extragalactic sources, which dominate 
above the ankle. This can be splitted into the flux from the GC source, dJ/dil{E), 
and the contribution from all the other sources, dJbkg / dil{E) , i.e. dJtot/dfl — 
dJbkg/dn + dJ/dfl. 

The total (measured) differential energy spectrum of CRs can be approximated 



in the range of energies from 4 x 10^^ eV < E < 6.3 x lO'^^ eV, and C = 
(9.23 ±0.65) X 10-33m-2s-isr~ieV~\ 

The intensity of the total flux received at the Earth from the source is J{E) — 
Jq{E)A{E), where A{E) is the amplification computed in the previous section and 
for the actual source flux we can consider a power law dependence with energy 
Jo{E) (X i?"^-^, which is typical of Fermi acceleration in shocks. 

If we assume that the main contribution to the anisotropy comes from the 
GC neutron source, and hence consider that the background CR flux is essentially 
isotropic, we can write the observed anisotropy as A(£') = As{E)f{E), where 
f{E) = -Fb(s)/^o(6fcg) is the ratio of the contributions to the right ascension monopole 
from the source flux and the background flux. For an anisotropic flux, as in this case, 
f{E) depends on the exposure of the experiment. The anisotropy of the source flux, 
As{E), can be computed from the harmonic analysis of the right ascension distribution 
of the GC events as discussed in Section 13.11 We show in Figure the expected 
anisotropy as a function of the energy for different values of the turbulent magnetic 
field amplitude and exposures. For a random field amplitude Brms = 3 ^G, the 
anisotropy grows sharply around E ~ lO'^^'* as observed by AGASA. This is not 
the case for smaller values of the rms amplitude of the turbulent field as is apparent 
from the curve labeled Brms — 1 fJ-G in the Figure, as the CR arrival directions are 
not spread sufficiently at small energies to wash out the anisotropics. Notice that 
the effects associated to a harder GC source spectrum are partially compensated by 



by CHI 
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Figure 5. Anisotropy as a function of energy for random magnetic fields of 
strengths Brma = 1 and 3 fiG, and for exposures corresponding to full sky (FS), 
AGASA (AG) and for a southern observatory at b = —35° (S). 



the effects due to the increased amphfication of the fluxes for decreasing energies, 
leading to a more moderate energy dependence. We have normahzed the source flux 
in such a way that the anisotropy is 4% at i? ~ lO^^eV for the AGASA exposure. 
This corresponds to a fraction / ~ 0.026 contributed by the GC source to the right 
ascension monopolc at that energy. 

The sharp faU in the anisotropy appearing at _E 2± lO^^ "^ eV for the AGASA 
curve arises because above that energy most CRs arrive as direct neutrons from the 
GC direction, which is out of the held of view of AGASA. The possible existence of a 
cutoff in the acceleration energy of the source around that range of values is actually 
not necessary to explain AGASA's results. 



3-4-. Source luminosity 

We can estimate the luminosity of the GC source which corresponds to the 
normalization of its flux at the Earth adopted in the previous Section to flt the 4% 
amplitude of the right ascension flrst harmonic measured by AGASA. This is given by 

Lac^^^rl J dEEJoiE), (7) 

where rg is the distance to the GC, and the source flux in the absence of magnetic 
fields is given by JoiE) = J{E)/A{E). The source flux J{E) at lO^^eV can be related 
to the background flux J^kg at that energy using the value of f[E) obtained in the 
last Section. With the amplification computed in Section rOl the local fiux at lO'^^eV 
results Jo(10^* eV) 1.3 x 10~'^°TO~^s~^eV~^. For a power law energy spectrum of 
the source Jo{E) oc E~'^ '^, the total luminosity in the decade of energy from 10^^'^ 
to lO^S-S then results LgcCIO^"^'^ - lO^^-^) ~ 4 x lO^^ erg/s. This value is more than 
an order of magnitude lower than the maximum luminosity (10^^ erg/s) estimated for 
Sgr A* (8j mentioned in the Introduction. 

The expected flux of direct neutrons is a factor e-''o/T=^" of the source flux Jq{E). 
The integral of this in the energy range from 10^^'^ eV to 10^'*'^ eV corresponds then 
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Figure 6. Ratio of tlie density in each direction to the mean density at that 
declination for E = 1 EeV and Brms = 3 ^G. 

to a direct neutron flux of 2 x 10~^'^m~^s~^, which is not very different from the value 
(9 ± 3) X 10~^^m~^s~-^ reported by SUGAR for their point-Uke excess 

In order to give an idea of the distribution of the excess of events over the ones 
expected for an isotropic distribution, we made a map (similar to those constructed by 
AGASA and SUGAR) in which we added to an isotropic background the CR density 
expected from a GC neutron source of the required intensity to fit the AGASA's right 
ascension anisotropy. The map in Figure displays the ratio of the intensity at each 
direction to the mean density at that declination (following AGASA's analysis). We 
have smoothed the map with a window of 3° radius to take into account the typical 
angular resolution at these energies. Although the particular details of the map will 
depend on the GMF structure and strength, that are poorly known, we think that 
the main characteristic of the distribution are picked with the realistic model adopted. 
Notice however that the proton signal around the GC is sensitive to the actual location 
of the reversals in direction of the regular magnetic field. In the model adopted, in 
which the nearest reversal to the GC direction is at less than 1 kpc from the Earth, 
the excess observed at southern galactic latitudes is more significant than the one at 
northern latitudes, but this could change if the reversal was at a smaller galactocentric 
radius. Also, a vertical component of the magnetic field could displace these excesses 
in the direction parallel to the galactic plane. 

4. Conclusions 

We examined the possibility that a source of neutrons at the GC could be responsible 
for the anisotropy measured by AGASA and SUGAR. As the GC itself is out of the 
field of view of AGASA, this observatory would not detect direct neutrons from a 
source in the GC, but it would just see the protons produced by the neutron decays. 
As their trajectories are bent by the GMF, the distribution of protons arriving to 
the Earth depends strongly on the energy of the particles and on the structure and 
strength of the field, which are not well known. We assumed a particular model for the 
regular field, which is consistent with the main observational evidences, and considered 
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different reasonable amplitudes for the turbulent magnetic field. We found that there 
is a sharp rise in the anisotropy at energies around 1 EeV, as observed by AGASA, 
for the largest amplitude of the turbulent magnetic field considered, Bj-ms = 3 ^G. At 
energies larger than 2 EeV, essentially no particles are expected from the GC source 
in the field of view of AGASA, what explains the disappearance of the anisotropy 
independently of the source acceleration cutoff. The predictions are however very 
different for an observatory in the south. This one should detect the strong point-like 
signal from the direct neutrons in addition to the extended signal from the protons. 
At energies below 1 EeV the neutron signal is suppressed due to the reduced decay 
length of the neutrons and the expected anisotropy has a similar behavior as that 
of a northern observatory. At larger energies, and up to the maximum acceleration 
energy of the source, the point-like signal should be clearly seen. It is interesting 
to notice that the neutron source luminosity needed to fit the 4% right ascension 
anisotropy detected by AGASA (due to the protons in this model) give rise to a 
point-like neutron flux consistent with the one measured by SUGAR. However, unless 
there is some systematic pointing error in their data, the offset of the signal from the 
GC direction cannot be explained. Moreover, for the GMF adopted, the phase of the 
first harmonic in right ascension obtained (~ 330°) is somewhat larger than AGASA's 
one. Another important point is that within this scenario the excess observed along 
the Cygnus region is also related to the central source (for this it is crucial that the 
regular magnetic field be along the spiral arm and not just azimuthal). This excess 
disappears for E > 2x 10^^ eV due to the impossibility for the GMF to confine the CRs 
along the spiral arms, independently of the source upper energy cutoff. The alternative 
explanation of a diffusive flux of charged particles from the GC direction [HlEIj should 
give rise to a wide angle anisotropy, with a deficit of particles in the galactic anticenter 
direction and a cosine angular dependence with respect to the direction of maximum 
flux, quite different from the neutron source signal. More data from an observatory in 
the south would be very helpful to confirm or reject these hypothesis. Although the 
AUGER observatory ^HI in construction in Argentina is focused to the study of CRs 
with energy above 10^^ eV, it may be able to detect lower energy events (down to 
10^^ eV). A dedicated detector for the 10^^ eV to 10^^ eV energy range is also under 
study 

If the sharp peak in the GC direction is seen, confirming the neutron source 
hypothesis, then the distribution of the protons would give information about the 
structure of the GMF. 
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